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摘 要 : 奶牛 步 态 时 相 是 反映 奶牛 健康 及 跨行 严重 程度 的 重要 指标 。 为 准确 自动 识别 奶牛 步 态 时 相 ， 本 研 


究 提 出 一 种 融合 高 


混合 模型 (Gaussian Mixture Model, GMM) 和 隐 马 尔 科 夫 模型 


Hidden Markov Model, 


HMM) 的 无 监督 学 习 奶 牛 步 态 时 相识 别 算法 GMM-HMM。 使 用 惯性 测量 单元 采集 奶牛 后 肢 加 速度 和 角速度 
信号 ， 通 过 卡尔 曼 滤 波 消除 噪声 ， 筛 选 并 提取 特征 值 ， 构 建 CGMM-HMM 模型， 实现 奶牛 静 立 相 、 连 续 步 态 


中 的 站 立 相 和 摆动 相等 3 种 步 态 时 相 的 自动 识别 。 结 昌 


89.28%. 90.95% 和 90.91% ， 连 续 步 态 中 的 站 立 相 识 昂 


89.06%， 连 续 步 态 中 的 摆动 相识 别 的 准确 率 、 召 回 率 和 分 别 为 86.67%、91.51% 和 89.03%。 奶 牛 步 态 分 割 


表明 ， 静 立 相识 别 的 准确 率 、 召 回 率 和 石 ,分 别 为 


的 准确 率 、 召 回 率 和 分 别 为 91.55%、86.71% 和 


的 准确 率 为 91.67%， 相 较 于 基于 事件 的 峰值 检测 法 和 动态 时 间 规 整 算法 准确 率 分 别提 高 了 4.23% 和 1.1%。 
本 研究 可 为 下 一 步 基于 穿戴 式 步 态 分 析 的 奶牛 跨行 特征 提取 提供 技术 参考 。 
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1 引言 


奶牛 踊 行 是 由 于 蹄 底 溃 疡 、 趾 间 皮 炎 和 腐 足 
病 ”等 各 种 肢 蹄 病 引起 奶牛 不 适 ， 导 致 其 承重 
能 力 不 足 以 及 步 态 异 常 。 奶 牛 踊 行 与 饲养 方式 、 
牧场 环境 卫生 等 有 关 ”， 影 响 奶牛 正常 生产 生 
活 ， 进 而 导致 产 奶 量 下 降 、 个 体 产 奶 年 限 缩短 以 
及 治疗 成 本 上 升 等 问题 ， 严 重 影响 牧场 经 济 效 
益 。 当 前 中 国 奶牛 肢 蹄 病 发 病 率 较 高 KITEE 
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过 了 30% ”， 人 工 识 别 存 在 效率 低 .主观 判断 存 
在 偏差 .劳动力 成 本 较 高 和 容易 造成 奶牛 应 激 反 
应 等 问题 ， 发 展 奶牛 踊 行 自动 识别 技术 需求 
强烈 “。 

奶牛 幢 行 步 态 分 析 的 关键 是 在 奶牛 连续 步 态 
中 识别 摆动 相 和 站 立 相 等 步 态 时 相 ， 通 过 分 析 奶 
牛 各 个 步 态 时 相 的 运动 强度 和 持续 时 间 判 断奶 牛 
由 行程 度 。 奶 牛 的 步 态 满 足 一 定 的 周期 性 和 规律 


基金 项 目 : 国家 自然 科学 基金 项 目 (32102600); 中 国 农 业 科学 院 科 技 创新 工程 项 目 (CAAS-ASTIP-2016-AI); 中 央 级 公益 性 科研 院 


所 基本 科研 业务 费 专 项 (JBYW-AII-2021-33, JBYW-AII-2022-36) 


作者 简介 : 张 楷 (1997 一 )， 男 ,硕士 研究 生 ， 人 研究 方向 为 冀 禽 行为 识别 与 畜牧 物 联网 技术 。E-mail: 82101212474@caas.cn 
“通信 作者 : 刘 继 芳 (1965 一 )， 男 ,博士 ,研究 员 ， 研 究 方向 为 农业 信息 管理 。E-mail: liujifang@caas.cn 


54 


ChinaXiv 合 作 期 刊 


202302.00167v1 


chinaXiv 


性 ， 但 在 现实 环境 中 ， 奶 牛 步 态 受 环境 和 个 体 差 
异 的 影响 较 大 ， 自 适应 奶牛 步 态 分 析 需 要 针对 现 
实情 况 采 取 相 应 的 技术 方法 。 现 较 多 采用 计算 机 
视觉 技术 和 可 穿戴 式 设备 中 来 进行 奶牛 步 态 分 
析 。 计 算 机 视觉 技术 具有 非 接触 、 免 应 激 、 低 成 
本 、 高 通 量 、 不 干扰 奶牛 正常 活动 等 优点 “, 已 
被 用 于 奶牛 步 态 分 制 " 、 步 态 特 征 提取 以 及 
跨行 检测 ”中 。 但 计算 机 视觉 容易 受到 环境 光 
线 、 复 杂 背 景 以 及 奶牛 相互 诞 挡 的 影响 "。 可 
穿戴 式 传感器 O 具有 不 限于 特定 测量 区 域 以 及 
不 受 光照 和 遮挡 影响 的 优点 ， 目 前 ， 基 于 可 穿戴 
式 传 感 吉 进行 奶牛 跨行 的 研究 热点 集中 在 奶牛 步 
态 分 割 和 特征 提取 方面 。 常 见 的 奶牛 步 态 分 割 方 
式 主要 有 基于 事件 的 峰值 检测 法 和 基于 时 域 的 动 
态 时 间 规 整 方法 两 种 。 峰 值 检 测 法 最 先 应 用 于 奶 
牛 步 态 分 割 和 识别 ， 但 其 容易 将 一 些 非 步 态 
的 峰值 识别 为 步 态 事件 ， 导 致 低 准确 率 和 高 召回 
率 。 动 态 时 间 规 整 法 是 指 通过 计算 两 个 序列 的 相 
似 性 进行 步 态 分 割 ， 设 定 国 值 以 及 建立 最 优 匹配 
模板 ， 自 动 化 程度 不 高 1。 

隐 马 尔 科 夫 模 型 (Hidden Markov Model, 
HMM) 是 马尔 科 夫 模型 发 展 而 来 的 机 器 学 习 算 
法 ,核心 思想 是 通过 人 研究 观测 序列 来 辨识 不 可 见 
状态 序列 。 近 年 来 ， 许 多 研究 在 人 的 步 态 识别 上 
使 用 HMM。 在 人 体 步 态 的 研究 中 分 析 发 现 ， 步 
态 信 号 符合 马尔 科 夫 链 的 时 序 变 化 ， 针 对 时 序数 
据 的 动态 过 程 具有 重 棒 性 以 及 能 够 处 理 牧场 条 件 
下 奶牛 步 态 数据 中 的 异 质 性 ， 可 以 进行 步 态 阶段 
识别 呈 。 不 同步 态 分 割 方法 对 比 显 示 ， 基 于 
HMM 的 分 割 方法 优 于 动态 时 间 规 整 模 板 匹配 法 
以 及 峰值 检测 法 所。 由 于 本 研究 解决 的 是 长 时 
间 连 续 序 列 的 奶牛 步 态 时 相识 别 问题 ， 观 测 数据 
量 较 大 ， 所 以 引入 高 斯 混合 模型 (Gaussian Mix- 
ture Model，GMM)。GMM 可 以 无 监督 地 对 观测 
数据 进行 聚 类 ，GMM 和 HMM 混合 算法 曾 用 于 
人 类 步 态 识别 ， 采 用 无 监督 的 学 习 方 式 ， 得 到 帕 
金森 患者 步 态 分 割 ， 准 确 率 为 96.2% " 。 HAT, 
基于 GMM FI HMM 融合 的 奶牛 步 态 时 相识 别 研 


究 尚 未 展开 。 

综 上 所 述 ， 现 有 技术 方法 初步 实现 了 奶牛 步 
态 分 割 ， 但 还 未 深入 到 奶牛 步 态 时 相识 别 ， 并 且 
依旧 存在 数据 标注 困难 、 准 确 率 均 未 达到 实际 应 
用 要 求 的 问题 。 为 进一步 改进 奶牛 步 态 时 相 上 自动 
识别 方法 以 及 提高 其 准确 率 ， 本 研究 提出 一 种 融 
合 GMM 和 HMM 的 无 监督 学 习 奶 牛 步 态 时 相识 
别 算法 GMM-HMM， 利 用 惯性 测量 单元 采集 奶 
牛 步 态 数据 ， 结 合 视频 标注 构建 奶牛 步 态 时 相识 
别 数据 集 ; 筛选 出 冠状 轴 角 速度 的 一 阶 差 分 作为 
奶牛 步 态 时 相识 别 的 特征 值 ， 减 少 奶牛 个 体 差异 
对 识别 效果 的 影响 ， 并 分 析 奶 牛 步 态 分 割 效 果 ， 
为 奶牛 步 态 分 析 提 供 一 种 新 的 技术 手段 。 


2 材料 与 方法 


2.1 试验 设备 与 场地 


试验 在 中 国 农业 机 械 化 科学 研究 院 北京 农机 
试验 站 奶牛 养殖 场 进行 。 试 验 对 象 为 5 头 随 机 选 
择 的 荷 斯坦 奶 牛 ， 利 用 MEMS 微型 惯性 测量 单 
元 传感器 (LPMS-B2 型 号 ) 采集 奶牛 后 肢 加 速 
度 和 角速度 数据 ， 采 样 频率 为 50 Hz， 陀 螺 仪 测 
量 范围 为 +2000 dps。 在 奶牛 活动 场 内 ， 选 定 一 
条 长 20m、 宽 1.5m 走 廊 作 为 测量 场地 ， 测 量 场 
地 示意 图 如 图 1 所 示 。 使 用 佳能 EOS 5D Mark II 
相机 同步 拍摄 奶牛 通过 走廊 的 视频 。 


Tea 
区 域 


图 1 测量 场地 示意 图 


Fig. 1 Schematic diagram of measuring site 


Bi 
up 


ChinaXiv 合 作 期 刊 


55 


202302.00167v1 


chinaXiv 


2.2 数据 采集 与 预 处 理 


2.2.1 数据 采集 

由 于 奶牛 蹄 病 发 生 于 后 肢 的 概率 较 高 一 ， 
以 及 为 方便 收集 数据 ， 统 一 将 传感器 绑 在 奶牛 右 
后 胶 上 。5 头 试验 奶牛 佩戴 方向 完全 相同 ， 在 传 
感 器 中 ， 定 义 x 轴 为 奶牛 前 进 方向 ，y 轴 为 竖 直 
方向 ; z 轴 为 冠状 轴 方 向 ， 如 图 2 所 示 。 奶 牛 在 
测量 场地 自由 活动 ， 相 机 以 50 fs 帧 率 录制 奶牛 
通过 走廊 的 视频 。 根 据 传 感 器 的 时 间 戳 与 相机 记 
录 的 视频 时 间 实 现 图 像 数 据 与 传感器 数据 时 间 同 
步 。 试 验 中 ， 采 用 人 工 引 导 的 方式 确保 所 有 试验 
奶牛 多 次 通过 走廊 ， 记 录 充 足 的 图 像 数据 。 将 人 
工 标注 的 视频 数据 作为 奶牛 步 态 时 相识 别 的 黄金 
标准 。 试 验 采 集 了 2022 年 2 月 23 日 14:00 一 16:30 的 
奶牛 步 态 数据 ， 每 个 传感器 有 240,000~260,000 
条 数据 记录 ， 拍 摄 的 奶牛 通过 走廊 的 视频 共计 
18 段 。 


a 


图 2 传感器 佩戴 方式 及 传感器 坐标 


Fig. 2 Sensor wearing method and sensor coordinates 


2.2.2 ”数据 预 处 理 

LPMS-B2 传感器 数据 中 与 奶牛 步 态 的 运动 
学 数据 直接 相关 的 是 加 速度 和 角速度 ， 使 用 3 轴 
加 速度 和 3 轴 角 速度 数据 为 基础 进行 奶牛 步 态 时 
相识 别 。 由 于 传感器 同时 测量 9 轴 数 据 ， 高 频 噪 
声 较 少 ， 噪 声 多 受 系 统 干扰 的 影响 ， 为 去 除 系统 
的 噪声 和 干扰 影响 ， 还 原 真 实 的 动态 数据 ， 采 用 
卡尔 曼 滤 波 (Kalman filtering) 进行 数据 预 处 
理 。 本 试验 使 用 的 传感器 集成 加 速度 计 、 陀 螺 仪 
和 磁力 计 三 个 部 分 ， 同 时 测量 x、y、z 轴 三 个 方 
回 的 线 加 速度 、 角 速度 以 及 磁场 强度 。 在 测量 过 
程 中 合计 9 轴 的 数据 互相 影响 并 且 单 位 不 同 ， 使 
用 卡尔 曼 滤 波 可 以 减轻 传感器 噪声 ， 降 低 系统 误 


差 和 测量 误差 。 该 传 感 名 内置 滤波 器 ， 在 采集 数 
据 之 前 已 将 滤波 器 设 定 完成 ， 得 到 的 数据 噪声 因 
素 已 经 消除 ， 可 以 直接 利用 数据 进行 模型 分 析 。 
将 图 像 数 据 采 用 FFmpeg 软件 进 行 视频 分 帧 处 理 ， 
通过 逐 帧 比 对 的 方式 ， 确 定 蹄 离 地 和 蹄 落地 两 个 
步 态 事 件 ， 再 与 传感器 数据 进行 匹配 。 由 图 2 可 
知 ，x 轴 的 加 速度 可 以 表示 奶牛 基本 的 运动 学 特 
征 ，z 轴 方向 为 水 平 且 垂直 于 前 进 方向 向 外 ，z 
轴 角 速度 可 以 捕 换 奶牛 细微 的 牛 蹄 运动 ， 可 以 进 
一 步 分 析 奶 牛 的 步 态 情况 ， 所 以 分 析 奶 牛 后 胶 的 
x 轴 加 速度 和 z 轴 角速度 可 以 反应 奶牛 步 态 事件 ， 
更 容易 实现 奶牛 步 态 事件 与 传 感 融 数据 的 同步 。 
x 轴 加 速度 和 z 轴 角速度 折线 图 对 应 的 奶牛 步 态 
事件 如 图 3 所 示 。 


2.3 训练 集 和 测试 集 数据 构建 


数据 集 构建 以 测试 集 为 基准 ， 测 试 集 构建 时 
使 用 人 工 标 注 的 方式 ， 将 图 像 数据 与 传感器 数据 
进行 匹配 ， 在 传 感 融 数据 中 标 出 步 态 事件 〈 蹄 落 
地 ， 蹄 离 地 )， 图 像 数据 是 随机 采集 的 ， 符 合 测 
试 集 构建 的 规则 。 观 察 6 轴 传感器 数据 并 做 了 大 
量 尝试 后 发 现 ， 需 要 找 出 周期 性 最 明显 的 数据 ， 
方便 数据 标注 也 更 容易 找到 特征 值 ， 最 终 确定 将 
z 轴 的 角速度 作为 本 次 试验 数据 ， 如 图 4 所 示 。 

为 消除 佩戴 传感器 时 奶牛 行为 以 及 刚 佩 戴 好 
传感器 后 奶牛 不 适应 带 来 的 试验 误差 ， 将 原始 数 
据 的 前 80,000 个 记录 点 ( 约 26.7 min) 的 数据 作 
为 误差 项 去 除 ， 同 理 ， 将 后 80,000 个 记录 点 同时 
去 除 。 由 于 传感器 存在 时 间 惟 ， 利 用 视频 时 间 可 
以 实现 奶牛 步 态 时 间 与 传感器 数据 同步 ， 每 头 牛 
有 80~-90 min 的 传感器 数据 ， 提 取 有 视频 标注 的 
传感器 数据 作为 测试 集 ， 剩 余 的 传感器 数据 去 除 
误差 项 作为 无 监督 的 训练 集 ， 其 中 包含 每 头 牛 至 
少 30 min 以 上 的 高 质量 传感器 数据 ， 测 试 集 中 连 
续 步 态 有 162 步 ， 训 练 集中 连续 步 态 有 729 步 。 


2.4 GMM-HMM 融合 算法 描述 


本 人 研究 采用 GMM-HMM 融合 算法 进行 奶牛 
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注 :AX 为 X 轴 加 速度 ,GZ 为 z 轴 角速度 


图 3 x 轴 加 速度 和 z 轴 角速度 折线 图 对 应 的 奶牛 步 态 事件 


Fig. 3 Gait events of dairy cows corresponding to the curve of x-axis acceleration and z-axis angular velocity 


步 态 时 相识 别 研究 。 

奶牛 步 态 时 相 为 奶牛 行走 的 单个 步 态 周期 中 
不 同 阶 段 ， 本 人 研究 将 连续 步 态 分 为 摆动 相 和 站 立 
相 。 摆 动 相 为 奶牛 蹄 离 地 到 蹄 着 地 的 过 程 ， 站 立 
相 为 蹄 着 地 到 蹄 离 地 的 过 程 。 在 连续 的 步 态 序列 
中 ， 摆 动 相 和 站 立 相 是 交替 出 现 的 ， 但 是 本 试验 
过 程 中 ， 奶 牛 除了 行走 以 外 ， 还 存在 静止 站 立 的 
状态 。 因 此 ， 本 研究 将 奶牛 四 胶 同 时 着 地 静止 不 
动 的 状态 定义 为 静 立 相 ， 区 别 于 连续 步 态 序列 
中 ， 奶 牛 四 肢 交 换 站 立 的 站 立 相 。 静 立 相 、 摆 动 
相 和 站 立 相 都 可 以 反映 奶牛 的 运动 状态 ， 同 时 奶 
牛 步 态 时 相 的 识别 可 以 实现 对 奶牛 的 步 态 分 割 ， 
识别 的 连续 站 立 相 和 摆动 相 就 是 一 个 完整 的 
步 态 。 

GMM 是 多 个 高 斯 概率 密度 函数 的 线性 组 合 。 
为 满足 奶牛 步 态 时 相识 别 的 鲁 棒 性 ， 需 要 奶牛 长 
时 间 佩 戴 传 感 涡 ， 采 集 的 测量 数据 包含 静 立 相 , 
和 连续 步 态 中 的 摆动 相 5, 和 站 立 相 5S,， 总 共 3 个 
步 态 时 相 ， 每 个 时 相对 应 一 个 高 斯 分 布 ， 测量 数 
据 符合 这 3 个 时 相 的 高 斯 混合 分 布 。 使 用 该 模型 
的 目的 是 通过 训练 找到 相应 参数 69， 进而 将 高 斯 


混合 分 布 分 复 ， 得 到 原始 的 每 一 个 高 斯 分 布 ， 实 
现 类 似 分 类 的 效果 ， 进 而 将 不 同 的 状态 分 类 。 高 
斯 混合 模型 的 实现 过 程 即 为 参数 0 的 估计 ， 采 用 
期 望 最 大 化 (Expectation Maximization, EM) 算 
法 进行 参数 估计 ""， 该 方法 基于 极 大 似 然 估计 
算法 可 表示 为 : 


P(X0)- >oe(X) (1) 
1 | 1 
@(X) = ——_—_,, -~(Xx- 
(xX) E exp) -z 
wj- (2) 
0= {8 05); i= 1,2,.%,n (3) 


其 中 , 2 是 GMM 的 阶 数 ， 本 研究 中 n=3; mw， 
是 权 值 ， 且 w,> 0, >w,= 1; B(X) 是 第 i 个 时 


相 的 高 斯 分 布 概率 密度 函数 ; D, u, SAP GE 
据 的 维 数 、 期 望 和 协 方差 。 

HMM ”由 马尔 科 夫 模型 发 展 而 来 ， 是 一 种 
动态 贝 叶 斯 的 生成 模型 以 及 有 向 图 模型 ， 由 隐 状 
态 序 列 和 观测 序列 构成 ， 其 中 每 个 观测 数据 代表 
一 个 状态 。 核 心思 想 是 由 于 状态 序列 的 不 可 直接 
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注 :AX、AY、AZ 分 别 为 X、y、zZ 轴 的 加 速度 ;GX、GY、GZ 分 别 为 Xx、y、z 轴 的 角速度 
图 4 传感器 6 轴 数 据 折 线 图 


Fig. 4 Sensor 6-axis data line chart 


观测 的 性 质 ， 通 过 研究 观测 序列 来 辨识 不 可 见 隐 
状态 序列 ， 如 图 5 所 示 。 
HMM 由 三 个 参数 构成 : 初始 状态 概率 x、 


tO E 
wm @) @-@ © ~ 
注 :5 为 隐 状 态 ,0O, 为 观测 状态 ,i=1,2,3………n 


图 5 隐 马 尔 科 夫 模型 结构 
Fig. 5 Hidden Markov Model structure 


状态 转移 概率 4 以 及 发 射 概率 〈 隐 状态 发 生 观 测 
事件 的 概率 ) B。z 和 4 构成 了 隐 状 态 之 间 的 状态 
转移 ， 如 图 6 所 示 。 其 中 ，5 为 静止 相 ，S,、S; 为 
连续 步 态 中 的 摆动 相 和 站 立 相 ，a; 为 第 i 个 状态 
转移 到 第 ;个 状态 的 概率 i，j=1，2，3。B 将 隐 状 
态 和 观测 序列 联系 ， 观 测序 列 由 状态 序列 生成 。 
该 参数 集 可 以 表示 为 : 

A = {n, A, B} (4) 

HMM 可 以 解决 三 个 问题 : 

(1) 评估 问题 。 给 定 观测 序列 和 参数 集 4， 
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图 6 隐 马 尔 科 夫 模型 状态 转换 的 过 程 


Fig. 6 State transition processof Hidden Markov Model 


计算 观测 序列 出 现 的 概率 ， 采 用 前 向 -后 向 (For- 
ward-Backword) 算法 实现 。 

(2) 学 习 问 题 。 给 定 观 测序 列 ， 佑 计 模 型 参 
数 集 4， 方 法 的 核心 是 使 得 观测 序列 出 现 的 概率 
最 大 的 参数 即 为 模型 参数 ， 采 用 极 大 似 然 估计 的 
方法 (EM 算法 或 者 Baum-Welth 算法 ) 实现 。 

(3) 解码 问题 。 已 知 参数 集 4 和 观测 序列 ， 
计算 概率 最 大 的 状态 序列 ， 采 用 Viterbi 算 法 实 
现 。Viterbi 算 法 本 质 上 是 一 个 动态 规划 的 最 短路 
问题 ， 找 到 使 得 状态 序列 发 射 概率 最 大 即 由 观测 
状态 序列 最 可 能 对 应 的 状态 序列 。 

在 本 研究 中 ， 先 进行 学 习 问 题 的 训练 ， 利 用 
训练 集 得 出 参数 集 ， 再 基于 得 到 的 模型 参数 ， 利 
用 测试 集 对 观测 序列 进行 解码 。 奶 牛 的 观测 数据 
为 观测 序列 ， 奶 牛 的 步 态 时 相 序列 为 状态 序列 ， 
步 态 时 相 为 隐 仿 状态。 为 解决 奶牛 步 态 时 相 分 割 
的 鲁 棒 性 ， 本 研究 针对 的 是 长 时 间 的 连续 观察 序 
列 ， 采 用 GMM 5 HMM 融合 模型 (GMM- 
HMM)， 此 时 参数 集中 的 发 射 概率 8 由 一 组 观察 
值 的 高 斯 分 布 密度 函数 蔡 代 ， 即 发 射 概率 8B 用 高 
斯 混合 模型 的 参数 进行 替代 ， 新 的 参数 集 为 : 

A = (a, A, 0, pt, X} (5) 

其 中 ，r、4 分 别 是 HMM 的 初始 概率 和 状态 
转移 概率 ; ww 和 5 分 别 为 观测 序列 高 斯 混合 分 
布 的 权重 ,均值 和 协 方差 。 

将 GMM-HMM 中 高 斯 分 布 数 设 定 为 3， 对 
应 静 立 相 、 扎 动 相 和 站 立 相 3 种 步 态 时 相 。 将 构 
建 好 的 数据 集 提取 特征 值 ， 将 z 轴 角速度 的 1 阶 
差分 作为 测试 集 带 和 模型， 进行 学 习 问 题 。 首 先 


对 参数 初始 概率 上 和 转换 概率 4 进行 初始 化 ， 然 
后 根据 EM 算法 反复 迭代 训练 得 出 参数 w、4 和 
区 ， 此 时 已 经 得 出 3 种 高 斯 分 布 ， 即 3 种 状态 ,x 
直接 按照 输入 的 训练 集 初 始 状态 定义 ， 再 根据 应 
用 统计 的 方法 计算 得 出 4， 公 式 如 下 : 
N; 

T 


i (6) 


其 中 ，a; 是 状态 i 转移 到 状态 j 的 概率 ; Ni、 
7 分 别 是 状态 i 转移 到 状态 j 的 个 数 、 各 个 状态 之 
间 转 换 的 总 个 数 。 

GMM-HMM 使 用 Python3.8 中 hmmlearn. 
hmm 模 块 实现 。 


2.5 奶牛 步 态 时 相识 别 


经 过 对 奶牛 行走 的 大 量 观察 以 及 查阅 资料 可 
知 ， 奶 牛 的 单个 步 态 周 期 为 1.2 s 左右。 为 满 
足 步 态 时 相识 别 的 自 适应 性 和 和 鲁 棒 性 ， 需 要 对 连 
续 长 时 间 序 列 的 步 态 数据 进行 研究 。 构 建 训练 集 
时 ， 考 虑 到 有 监督 学 习 的 训练 方式 需要 人 工 标注 
的 数据 量 过 于 庞大 且 效 率 低 ， 因 此 本 人 研究 采用 无 
监督 学 习 的 训练 方式 。 通 过 GMM-HMM 融合 算 
法 ， 从 传感器 的 运动 学 数据 中 找到 满足 条 件 的 特 
征 值 ， 对 连续 长 时 间 序 列 数据 进行 无 监督 的 聚 类 
以 及 参数 集训 练 ， 得 到 模型 参数 ， 用 测试 集 验证 
结果 ， 与 标注 数据 对 比 ， 得 出 识别 结果 ， 确 定 奶 
牛 步 态 时 相 ， 并 实现 奶牛 步 态 分 制 。 具 体 流程 如 
图 7 所 示 。 


a. = 


3 结果 与 分 析 
3.1 特征 值 选择 


对 奶牛 步 态 的 运动 学 数据 进行 统计 化 处 理 。 
将 2.3 中 选取 的 z 轴 角速度 采用 构建 滑动 时 间 窗 
和 不 使 用 构建 滑动 时 间 窗 的 方式 ， 提 取 数 据 的 最 
值 、 均 值 、 方 差 、 标 准 差 、1 阶 差 分 、 线 性 回归 
的 斜率 和 截 距 ， 以 及 平均 绝对 误差 。 最终 发 现 不 
加 窗 的 情况 下 ，1 阶 差分 可 以 很 好 地 区 别 奶 牛 步 
态 不 同时 相 的 数学 特征 。 图 8 为 不 同 特征 值 的 分 
布 比较 情况 ， 截 取 的 是 连续 步 态 周期 的 一 段 ， 两 
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EMA 


对 训练 集 进行 无 监督 学 习 的 聚 类 


奶牛 步 态 时 相识 别 


0 测试 集 


奶牛 步 态 分 割 


模型 参数 的 训练 


图 7 奶牛 步 态 时 相 的 自 适 应 识别 方法 流程 图 


Fig. 7 Flow chart of adaptive recognition method for gait phase of dairy cows 


侧 离 散 程 度 大 的 是 摆动 相 ， 两 条 黑 线 之 间 离 散 程 
度 小 的 是 站 立 相 。 由 图 8 可 知 ， 偏 度 在 站 立 相 和 


精度 不 足 ， 只 有 1 阶 差分 可 以 直观 且 高 精度 的 区 
分 站 立 相 和 摆动 相 ， 其 中 站 立 相 数据 较为 平缓 ， 


摆动 相 没 有 区 分 度 ， 最 值 和 均值 的 重合 度 过 高 且 摆动 相 数 据 波动 较 大 。 
没有 区 分 两 个 相位 的 直观 特征 ， 标 准 差 区 分 度 和 
。 “1 阶 差分 ' 
75 均值 
标准 差 
50 io - RA 5 
B 偏 度 for + 
Z 25 a g ae 
Bas ; a 区 A o ers oa 
0 Risa LT Tt hearst to ax 
-25 a . : : . 
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+ grep is 
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图 8 奶牛 步 态 不 同时 相 的 数学 特征 


Fig. 8 Mathematical characteristics of different gait phases of dairy cows 


#1 GMM-HMM 
Table 1 GMM-HMM M 


3.2 模型 参数 训练 结果 


模型 参数 


odel parameters 


经 过 训练 得 到 的 模型 参数 包括 不 同时 相 的 初 


S 时 相 初始 概率 权重 均值 方差 
始 概率 、 高 斯 混合 模型 的 权重 、 均 值 和 方差 ， 如 menm a i Too mem 
表 1 所 示 ， 公 式 (7) 表示 步 态 时 相 的 转移 概率 站 立 相 0 1 -0.251 3.904 
乍 阵 4。 转 移 概 率 和 矩阵 按照 静 立 相 、 站 立 相 、 摆 摆动 相 1 1 0.357 1976.724 


动 相 的 顺序 排列 ， 表 示 时 相 之 间 相 互 转换 的 
概率 。 
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0.987 0.012 0 
A=!0.064 0.904 0.031 (7) 
0 0.046 0.953 


从 所 得 参数 的 结果 可 知 ， 初 始 概率 即 为 训练 
集 的 初始 时 相 的 概率 ， 模 型 构建 时 以 摆动 相 为 初 
始 输入 时 相 ， 满 足 随 机 性 。 三 个 时 相 的 高 斯 混合 
模型 的 权重 都 为 1， 说 明 每 一 阶段 单独 表示 相应 
时 相 ， 没 有 时 相 之 间 的 交叉 影响 ， 提 高 了 模型 的 
可 解释 性 和 准确 性 。 转 移 概率 的 数值 也 表明 了 模 
型 构建 的 准确 性 。 每 一 组 均值 和 方差 代表 一 个 
时 相 。 


3.3 奶牛 步 态 时 相识 别 结果 


参数 训练 完成 后 进入 解码 阶段 。 将 测试 集 带 
入 GMM-HMM 识别 奶牛 步 态 时 相 ， 得 到 的 部 分 
测试 结果 如 图 9 所 示 。 图 中 红色 、 绿 色 和 蓝 色 部 
分 为 奶牛 的 静 立 相 、 站 立 相 和 摆动 相 。 从 图 9 可 
以 直观 发 现 ， 较 为 平缓 的 是 奶牛 没有 行走 的 静 立 
相 ; 站 立 相 是 奶牛 蹄 着 地 到 蹄 离 地 的 这 一 过 程 ， 
其 角速度 与 图 中 体现 的 一 致 ， 从 蹄 着 地 开始 减少 
再 增加 直到 蹄 离 地 ;连续 步 态 的 站 立 相 和 摆动 相 
交 蔡 出 现 ， 摆 动 相 的 角速度 增加 到 一 个 步 态 周期 
的 最 大 值 ， 然 后 减少 ， 直 到 蹄 落地 开始 受 力 
为 止 。 


,oo| 一 一 静 立 相 
一 一 站 立 相 
150| 一 一 摆动 相 


速度 //(deg:s !) 
3 


0.0 2.5 5.0 7.5 


10.0 12.5 15.0 17.5 20.0 
时 间 /s 


图 9 奶牛 步 态 时 相识 别 结果 


Fig. 9 Gait phase recognition results of dairy cows 


采用 精确 率 (Precision, P), 4 E% (Re- 
call, R), {$ X (Specificity, Sp), UR PE 
(Sensitivity, Se) 和 局 值 为 评价 指标 。 本 方法 的 
步 态 时 相识 别 结 果 中 每 一 条 测试 集 数据 都 会 得 到 
一 个 状态 ， 根 据 测 试 视 频 标 注 结果 来 验证 测试 
集 ， 结 果 如 表 2 所 示 。 

由 表 2 结 果 可 知 ，3 种 时 相 在 本 研究 方法 的 
已 值 相 差 不 大 ， 静 立 相 的 已 值 较 高 为 90.91%， 
识别 效果 最 好 ,说 明 静 立 相 与 连续 步 态 的 两 个 时 
相 区 分 效果 较 好 ; 摆动 相 的 召回 率 和 敏感 性 最 
高 ， 都 为 91.51%， 说 明 摆 动 相 的 特征 较为 明显 ， 
识别 成 功率 最 高 ; 站 立 相 的 准确 率 和 特异 性 最 


表 2 GMM-HMM 的 奶牛 不 同步 态 时 相识 别 结果 
Table 2 GMM-HMM recognition results of different gait 


phases in dairy cows 


时 相 PI% RI% Sp/% Sel% F, /% 


静 立 相 89.28 90.95 93.83 90.95 90.91 
站 立 相 91.55 86.71 95.51 86.71 89.06 
摆动 相 86.67 91.51 93.58 91.51 89.02 


高 ， 分 别 为 91.55% 和 95.51%。3 个 相位 识别 结 
果 的 误差 来 源 主要 在 于 数据 标注 时 数据 质量 以 及 
试验 环境 因素 ,例如 ， 视 频 清晰 度 不 足以 准确 判 
断奶 牛 蹄 离 地 、 蹄 落地 这 些 事件 引起 的 误差 以 及 
搜集 数据 时 牛 场 湿 软 地 面 引起 奶牛 二 次 发 力 以 及 
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站 立 不 稳 等 因素 。 

通过 检索 文献 发 现 ， 大 部 分 针对 奶牛 步 态 的 
研究 关注 奶牛 步 态 分 割 方法 为 提取 跨行 奶牛 步 态 
时 相 特 征 ， 本 研究 更 进一步 研究 奶牛 步 态 时 相识 
别 算法 。 试 验 通 过 人 工 引 导 奶 牛 的 方式 ,采集 了 
奶牛 在 不 同行 走 速度 的 步 态 惯性 测量 数据 ， 结 
表明 本 研究 提出 的 方法 能 够 克服 行走 速度 对 步 态 
时 相识 别 结果 的 影响 ， 准 确 识别 步 态 时 相 ， 有 望 
实现 非 约 束 环境 下 奶牛 步 态 时 相 的 自 适 应 识别 。 


3.4 奶牛 步 态 分 割 结果 对 比分 析 


针对 奶牛 连续 步 态 阶段 中 的 站 立 相 和 摆动 相 
进行 奶牛 步 态 分 割 。 在 奶牛 连续 步 态 中 ， 站 立 相 
和 摆动 相 构 成 一 个 完整 的 步 态 周期 ， 根 据 连 续 步 
态 中 的 站 立 相 和 摆动 相交 替 出 现 ， 奶 牛 步 态 时 相 
识别 连续 步 态 阶段 中 的 站 立 相 和 摆动 相 的 结果 即 
为 奶牛 步 态 分 割 的 结果 。 奶 牛 步 态 分 割 时 不 再 以 
每 条 数据 为 标准 ， 使 用 成 功 分 割 的 步 态 周 期 为 单 
位 进行 结果 评价 ， 此 时 的 标准 为 人 工 标注 的 完整 
步 态 周期 。 将 成 功 识别 的 步 态 周期 设 定 为 正 样 
本 ， 连 续 步 态 中 验证 162 步 ， 设 定 误差 为 5 帧 
(+0.1s) 。 为 区 别 3.3 中 奶牛 步 态 时 相识 别 的 结 
果 ， 重 新 定义 TP* 为 准确 识别 单个 步 态 周期 的 个 
数 ，FP* 为 非 单个 步 态 周期 ， 但 被 识别 为 步 态 周 
期 的 个 数 ，FN* 为 未 被 成 功 识 别 的 步 态 周期 ， 其 
中 7TP* 为 121 个 ，FN* 为 35 个 ，FP* 为 11 个 ， 此 
时 TN* 无 意义 ， 得 到 准确 率 为 91.67%。 

同时 ， 将 本 研究 利用 传感器 的 奶牛 步 态 分 割 
方式 准确 率 结果 与 基于 图 像 处理 的 奶牛 步 态 分 割 
算法 ， 包 括 支持 向 量 机 、K 最 邻近 法 以 及 随机 森 
林 等 进行 对 比 。 峰 值 检测 法 奶牛 步 态 分 割 准确 率 
为 87.44%， 动 态 时 间 规 整 奶 牛 步 态 分 割 准确 率 
为 90.57% 5 ， 本 研究 奶牛 步 态 分 割 准 确 率 为 
91.67%， 不 仅 识 别 了 奶牛 步 态 时 相 ， 而 且 相 较 于 
峰值 检测 法 和 动态 时 间 规 整 分 别 将 奶牛 步 态 分 割 
准确 提高 4.23% 和 1.1%。 与 基于 图 像 处 理 的 奶牛 
步 态 分 割 算法 支持 向 量 机 、K 最 近邻 和 随机 森林 
的 准确 率 分 别 为 87.60%、83.90% 和 82.56%， 本 


研究 方法 分 别提 高 了 4.07%、7.7% 和 9.11%。 

本 试验 步 态 分 割 成 功 识 别 160 步 ， 蹄 落地 的 
总 误差 为 -14.9 s, 平均 每 步 误 差 0.093 s， 蹄 离 地 
的 总 误差 为 -7.2s， 平 均 每 步 误 差 -0.045 s, WE 
0.1s 的 误差 。 其 中 “-“ 表 示 滞 后 (视频 和 传 感 
器 都 为 50 Hz)。 设 立 的 正 样本 是 站 立 相 ， 站 立 相 
是 奶牛 蹄 着 地 到 奶牛 蹄 离 地 这 一 过 程 ， 蹄 着 地 时 
奶牛 会 出 现 站 立 不 稳 、 滑 行 ， 以 及 二 次 用 力 重 新 
站 稳 的 现象 ， 导 致 蹄 着 地 这 一 过 程 误差 较 大 (B 
站 立 相 开始 阶段 ) ， 模 型 对 蹄 着 地 这 一 事件 判别 
不 准确 ， 蹄 落地 后 ， 传 感 器 会 再 次 收 到 抖动 的 信 
号 ， 直 到 稳定 状态 ， 才 会 识别 站 立 相 ， 进 而 导致 
蹄 着 地 这 一 事件 识别 会 清 后 。 虽 然 蹄 离 地 这 一 事 
件 ， 不 会 有 二 次 受 力 的 情况 ,误差 会 大 幅 减 小 ， 
但 是 松软 、 烙 连 的 地 面 情况 ， 也 会 出 现 较 小 的 沛 
后 情况 。 


4 结论 


本 研究 利用 可 穿戴 式 惯性 测量 单元 采集 奶牛 
后 肢 的 3 轴 加 速度 和 3 轴 角 速度 数据 ,经 过 数据 
预 处 理 和 特征 值 杜 别 ,选择 z 轴 角速度 的 1 阶 差 分 
作为 特征 值 ,通过 CMM-HMM 融合 算法 进行 奶牛 
步 态 时 相 的 识别 。 结 果 为 奶牛 步 态 时 相 和 奶牛 步 
态 分 割 提供 了 理论 依据 ， 有 具体 结论 如 下 。 

(1) 提出 了 一 种 GMM-HMM fit 4 AY FCA E 
学 习 奶 牛 步 态 时 相识 别 算法 ,识别 出 了 奶牛 的 静 
立 相 、 连 续 步 态 阶段 的 摆动 相 和 站 立 相 ， 其 中 静 
立 相 识别 的 准确 率 、 召 回 率 和 五 分 别 为 89.28%、 
90.95% 和 90.91%， 连 续 步 态 中 的 站 立 相 识别 的 
准确 率 、 召 回 率 和 分 别 为 91.55%、86.71% 和 
89.06%， 连 续 步 态 中 的 摆动 相识 别 的 准确 率 、 召 
回 率 和 五 分 别 为 86.67%、91.51% 和 89.03%。 进 
一 步 识 别 了 奶牛 步 态 时 相 ， 且 识别 效果 较 好 ， 为 
后 期 奶牛 跨行 特征 提取 提供 了 技术 方法 。 

(2) 根据 连续 步 态 的 站 立 相 和 摆动 相 实 现 了 
奶牛 步 态 分 割 ， 准 确 率 达到 了 91.67%， 相 较 于 
基于 事件 的 峰值 检测 法 和 动态 时 间 规 整 算法 准确 
率 分 别提 高 了 4.23% 和 1.1%， 分 割 效果 较 好 。 

为 满足 非 约 束 条 件 下 步 态 识别 的 应 用 要 求 ， 
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Abstract: The gait phase of dairy cows is an important indicator to reflect the severity of lameness. IThe accuracy of available 
gait segmentation methods was not enough for lameness detection. In this study, a gait phase recognition method based on 
Gaussian mixture model (GMM) and hidden Markov model (HMM) was proposed and tested. Firstly, wearable inertial sensors 
LPMS-B2 were used to collect the acceleration and angular velocity signals of cow hind limbs. In order to remove the noise of 
the system and restore the real dynamic data, Kalman filter was used for data preprocessing. The first-order difference of the an- 
gular velocity of the coronal axis was selected as the eigenvalue. Secondly, to analyze the long-term continuous recorded gait se- 
quences of dairy cows, the processed data was clustered by GMM in the unsupervised way. The clustering results were taken as 
the input of the HMM, and the gait phase recognition of dairy cows was realized by decoding the observed data. Finally, the 
cow gait was segmented into 3 phases, including the stationary phase, standing phase and swing phase. At the same time, gait 
segmentation was achieved according to the standing phase and swing phase. The accuracy, recall rate and F, of the stationary 
phase were 89.28%, 90.95% and 90.91%, respectively. The accuracy, recall rate and F, of the standing phase recognition in con- 
tinuous gait were 91.55%, 86.71% and 89.06%, respectively. The accuracy, recall rate and F, of the swing phase recognition in 
continuous gait were 86.67%, 91.51% and 89.03%, respectively. The accuracy of cow gait segmentation was 91.67%, which 
was 4.23% and 1.1 % higher than that of the event-based peak detection method and dynamic time warping algorithm, respec- 
tively. The experimental results showed that the proposed method could overcome the influence of the cow's walking speed on 
gait phase recognition results, and recognize the gait phase accurately. This experiment provides a new method for the adaptive 
recognition of the cow gait phase in unconstrained environments. The degree of lameness of dairy cows can be judged by the 
gait features. 


Key words: dairy cow lameness; gait phase; gait segmentation; gaussian mixture model; hidden Markov model; Kalman filtering 
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